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AEROELASTIC PROBLEMS OP AIRPLANE DESIGN* 
By H. G. Kussner 


1. INTRODUCTION 


For thirty years aeroelastic problems in airplane design were 
practically unknown, owing to the low speed and rigid design method of 
aircraft at that time. The first contributions to the unsteady lifting 
surface theory appeared then, but they created only academic interest. 

In the meantime, the speed of the airplanes continued to increase, while 
the aerodynamic refinement dictated consistently thinner and slenderer 
wing configurations. As a result, the stiffness and the oscillation 
phenomena connected with it have proved to be a limit of the technique 
just like the strength. 

The aeroelastic processes are complicated because they encompass 
the whole airplane with its many parts and parameters. Aeroelastic ity 
is in a stage of rapid development of its theoretical and experimental 
methods. Its scope expands continuously. The development started with 
studies on wing and aileron flutter. They were followed later by studies 
on tail and tab flutter. At very high speed, two new phenomena appeared, 
fluttering of the skin and vibrations on approaching sonic velocity, which 
are based on the instability of compression shocks at the curved sur- 
faces. An extreme case of flutter is the static instability at zero fre- 
quency. Relevant also is the distortion of the wing by the aerodynamic 
forces and the buckling and distortion of the skin. The reversal of the 
aileron effect due to the wing distortion, while concerning only the 
control, can be treated with the same formulas . 

As regards the gust stress and the stability of flight motion, the 
stiffness of the structural components also becomes so much more signi- 
ficant as the wing shapes become thinner and slenderer. This is why 
these zones are now included in aeroelastic ity and are treated in part 
by the methods developed for flutter. However, their discussion would 
go beyond the confines of the present report. It also applies to the 
aeroelastic problems of the rotating wings, the supersonic propellers, 
and the fan and turbine blades . The literature on aeroelastic ity since 
19^5 is too voluminous to be enumerated here. Some more recent compre- 
hensive descriptions on airplane flutter are presented in the reports 
by Kussner (ref. 1), Broadbent (ref. 2), Garrick (ref. 3), and 

"Aeroelastische Erobleme des Flugzeugbaus, " Zeitschrift fur 
Flugwissenschaften, 3. Jahrgang, Heft 1, Jan. 1955, PP- 1-18. 
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van de Vooren (ref . 4 ) , and in the introductory article by Scanlan and 
Rosenbaum (ref. 5) , which also contains further references. 


2. INDICATIONS OF. AIRPLANE VIBRATIONS (FLUTTER) 


Flutter rests on the fact that the airplane as a vibratory system 
absorbs energy from the air stream. This is so much more easily possible 
as the Bpeed of the airplane is greater. Both the damping energy and 
the consumption energy from the airstream increase with the flutter 
amplitude. Depending upon which of these energies grows faster, it may 
result in mild flutter with constant amplitude or vicious flutter with 
amplitude increasing steadily to failure. 

Simplest of all are the natural vibrations of elastomechanical 
systems in vacuum. If the system, exists in an incompressible fluid with- 
out moving in it, the Kelvin impulses are additive; they are the mass 
forces of the covibrating fluid masses. In the atmosphere, the Kelvin 
impulses generally amount to only a few percent of the mass forces of the 
airplane, hence do not change its vibratory behavior appreciably. For 
these so-called static vibrations, all points of the system pass almost 
simultaneously through zero position. The nodal lines, that is, the loci 
of disappearing- vibration amplitude, are then so much more sharply pro- 
nounced as the damping energy is less. It is customary, although not 
exactly correct, to identify all static oscillation modes as "oscillations 
with one degree of freedom. " 

The static oscillation modes are not capable, as a rule, of taking 
energy from the airstream; they rather give off energy on' the airstream, 
hence are damped additionally. Absorption of energy is generally 
possible only when a static oscillation mode is modified by the air forces 
which are additionally created on the oscillating wing during forward 
motion. In that event, there Is no longer a sharp nodal line; the new 
oscillation has several degrees of freedom. However, it is possible that 
certain static oscillation modes with one degree of freedom take up energy 
from the airstream in certain speed ranges, hence flutter, especially at 
higher Mach numbers, and In the supersonic range. ' (Cf . Runyan (ref. 6), 
Cunningham (ref. 7) , Watkins (ref. 8), Weber and Ruppel (ref. 8a).) 

To maintain a static oscillation mode with constant amplitude calls 
for considerable energy input, even for an oscillating system in vacuum. 

The cause for it is the solid friction :of the many individual components 
on each other, of which the airplane is riveted together, welded, or 
bonded, as well as the friction of the control bearings and control cables. 
Solid test rods have about 100 times lower damping than airplane wings of 
the conventional type, but even in these the damp'ing is probably produced 
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by solid friction on microscopic cracks, grain bound ar ies, and such. 

Thus surface effects in the widest sense are involved. It is not likely 
that these damping losses can be treated differentially; they can be 
taken into account only integrally and in time-average for the whole 
airplane. Since it is largely a solid friction that is Involved, the 
damping energy per cycle at forced vibrations of a certain mode and 
amplitude in a wide range is. nearly independent of the oscillation fre- 
quency. According to the Navier -Stokes theorem, fluid friction would 
yield damping proportional to frequency, hence must be excluded. Struc- 
tural damping losses can be reliably determined only by measurements on 
the finished airplane. It is customary to allow for the damping energy 
in the theoretical flutter formula by introducing a constant damping 
phase angle gj_ in all elastic stresses of a certain oscillation mode i; 

for conventional airplanes, gj_ ranges from 0.01 to 0.10 (radians). 

In still air, a linearly moving wing can exceed its critical speed 
without fluttering; it is then in unstable equilibrium. Any shock 
against such a wing initiates the entire natural frequency spectrum, 
that is, with a certain distribution of amplitudes and phase angles. 

But it does not signify that this distribution is suitable for absorbing 
energy from the air stream. When many shocks act on different parts of 
the airplane, however, it is increasingly possible that accidentally the 
right distribution of the natural oscillation modes for absorbing energy 
from the air stream and initial flutter is excited. This has been observed 
in wind-tunnel tests with dynamically similar models. Many cases of flut- 
ter were observed only in flight through gusty air, whereby a great variety 
of different shocks are exerted on the airplane. From these observations, 
it follows that a certain minimum size of shock Is necessary to initiate 
flutter. The extent to which this rests on boundary-layer effects or 
overcoming solid friction is an open question. 

When the flutter of an airplane or of a model wing is recorded, the 
record often reveals considerable departures from the harmonic or 
"sinusoidal" vibration, which may be due to aerodynamic causes and to 
the aforementioned solid friction. For example, there are complicated 
control surface forms with internal balance, for which the pressures 
are not linearly dependent on the control angle. Such controls may even 
flutter in two different partially foreign frequencies simultaneously. 

If the angle of attack of a wing at landing or pullout becomes so great 
that the flow starts to separate, there Is no further clear coordination 
between angle of attack and pressure distribution in steady flow. A 
rotary motion of the wing with one degree of freedom is then able to take 
up energy from the airstream. A periodic state of oscillation is sus- 
tained, which, as a result of the complicated processes in the boundary 
layer, is not linear at breakway. This type of flutter has been known 
for a long time and has lately been investigated by Halfman, Johnson, 
and Haley (ref. 9 )* 
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The airplane engineer does not wish, the fully developed flutter) 
he would rather prevent it, though it may please some experimenters. He 
wants to prove that all physically possible flutter phenomena, and even 
those with very low amplitude, lie outside the permissible speed of his 
airplane design. For this purpose, the harmonic, linearized oscillation 
theorem with simplified assumptions regarding air forces and structural 
damping has been utilized successfully up to now) it is particularly 
suitable for stability studies with arbitrary small amplitude . The first 
attempt involved the problem of making this simplest case of aeroelastieity 
amenable to mathematical treatment and theoretical solution. This problem 
is already very difficult and only approximately solvable, as will be 
shown in the following. 


' 3. THEORY OF STATIC OSCULATIONS 


The airplane is replaced by an elastomechanical system of n point 
masses and their connecting elastic members, which satisfy the classical 
theory of elasticity. How this substitution is done win be disregarded 
here. The deformations of such a replacement system under given small 
forces and particularly the natural vibrations under the mass forces can 
be computed by assuming infinitely small displacements, hence lineariza- 
tion. A suitable aid for solving this linearized problem is In matrix 
form. Its application to real (self -ad Joined) eigen -value problems in 
physics and engineering is generally known. Such a problem is involved 
here. 


A one -row matrix corresponds exactly to a vector „P^, a square 
matrix to a tensor a iic , whereby the first subscript (i) is always 

coordinated to the rows, the second subscript (k) to the columns. In 
the customary matrix calculus the subscripts are omitted. All these 
subscriptless "direct" calculi, to which the vector and tensor analysis 
themselves belong, are too little amenable to expansion) another drawback 
is that the significance of many operation signs must be marked and the 
sequence of the factors must not he changed. For the representation of 
the deformations of the lifting surface by polynomials or surface har- 
monics as well as their numerical integration by means df these functions, 
the subscripts are indispensable. Therefore we shall use subscripts, 
apart from a temporary application of the vector analysis in the flow 
theory. With respect to factors with identical subscripts or exponents 
i, k, and 2, m is added up from 1 to n, so far as. the subscripts 
are not bracketed. This is the only rule that needs to be followed. 

Consider a plate-like wing of .little thickness divided into n tables 
or "torsion boxes" (cf. Williams and Mech (ref. 10))) are the forces 

and the normal displacements of the n centers of gravity i = 1 

to n of these tables. The matrix equations read then .. 
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P i “ a ik^k 

a ik 

= a ki 

(1) 

^k = b kZ P Z 

b kl 

= b Zk 

(2) 

where ^ one -row matrices with 

n 

terms, 

a.^ and h square 


matrices with, n2 terms j is the stiffness matrix and the 

deformation matrix; both are symmetric. By equations (l) and ( 2 ) 


a ik b kZ “ & ll 

0 when i ^ l 

1 when 1 = 1 


( 3 ) 


The deformation matrix b-^ is obtained by solving the linear 
equation system (l) with respect to £jj., when a 1v is known. For 
abbr e viat ion 


b 


ki 




00 


For harmonic oscillation, the displacement of* point k is 

= A k exp jvt J - (5) 


where denotes the complex oscillation amplitude and v the cyclic 

frequency. The mass force is therefore 


p k = M (k) 


* 5 * 

at 2 


= - M (k) A k v2 


exp jvt 


(6) 


with denoting the masses concentrated in point k. The mann er by 

which a given continuous mass distribution may be best replaced by n 
separate masses in n given points, will be shown later. If 



6 


MCA OM 1402 


equilibria is to exist, the sum of the forces in every point must 
disappear. Multiplying equation (6) by b^/v^ gives, in the absence 

of air forces, according to equations (2), ( 5 ), (6), the principal 
equation 3 

M i “ (tM) iA = 0 * = iA 2 - (7) 


where 


<»>ik - b u“(lO (8) 

is termed dynamic matrix and A = l/v^ real eigen value. Equation ( 7 ) 
contains an nth order algebraic equation for the n eigen values Aj_ 

of the elastomechanic replacement system, the so-called secular equa- 
tion. It is obtained by equating to zero the denominator determinant 
of the linear equation system ( 7 ). The secular equation is solved 
either by using the Graeff process or, if n is great, ^by iteration. 
First estimate an amplitude distribution of the fundamental fre- 
quency, then resolve the matrix equation 


VA l = (bM) lk A£ „ (9) 


with respect to A'. With this approximate value, compute a new ampli- 
tude distribution A” and repeat the process until the values of A 

and Aj_|A n do no longer vary. This gives the exact solution of the 
fundamental frequency mode A^^^ and of the eigen value A-[_. Now the 
orthogonality condition _ — 

M i A i A i^ = 0 (10) 


is valid for all harmonic oscillations. 

From equations ( 7 ) and (10) follows then, by the same iteration 
process, the first harmonic oscillations of the system and so forth, 
until all n eigen values A^ and oscillation modes A^W of the 

system are obtained. 
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The matrix method is particularly suitable with automatic electronic 
computers; they can compute all the aforementioned matrix operations and 

iterations very quickly. Even for n = = 55 points the operations 


take only a few hours. This assures a wide range of adaptation of the 
replacement system to the actual airplane. The assumption of elastic 
beams is no longer necessary; even plate-like systems can be treated. 
As a rule, it is more convenient to fix the stiffness matrix a^ 


according to the design data, as indicated by Williams (ref. 10). The 
elasticity theory 'of thin plates has been described by Reissner and 
Stein (ref. 11), Fung (ref. 12), and Theodorides (ref. 15). 


k. THEORY OF AIR FORCES OF THE OSCILLATING LIFTING SURFACE 


To resolve aeroelastic problems, a theory of unsteady surface is 
necessary. Some new advances in this direction are described in the 
following based on references 1, lb- , 15, and 1 6 . Assuming perfect, 
compressible fluid, the Euler equations for the velocity potential $ 
and pressure p read 


d<fc 

p 0dt +P=P 0 


p o c o 


2 A^> + = 0 

dt 


(H) 

( 12 ) 


P 0 is air density, c Q velocity of sound in still air. Equation ( 12 ) 

is applicable to very small pressure changes only. The velocity of the 
fluid is u = Ai> . The substantial derivation with respect to time is 


jL = 

dt St 


+ UA 


(13) 


Let 


u = v + u ' 


where v is the velocity of undisturbed flow along the positive x axi6 
and u * the interference flow. Because the flight speed v changes 
rather slowly, as a rule, v = constant can be assumed. 'Then 
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du _ dp' 

at “ dt 


Ao =■ Ad 1 


( 14 ) 


In aeroelastic problems it is sufficient to know the corresponding 
very small air forces, since the aim. is limited to exploration of the 
start of the flutter mechanism in order to avoid flutter*. Hence, the 
linearization by assuming jn'l « vj the substantial derivation is then 


d_ 

dt 





A. 

3x 


( 15 ) 


Equations (l4) and (15) must be inserted in equations (ll) and (12) . 
These equations and their solutions can be simplified by the following 
substitutions - 



CD = 


.IvL 


v(l - 0 2 ) 


K = 


vL 


c 0 (l - W) 


(16) 


X = Lx 



( 17 ) 


X, Y, Z are ordinary Cartesian coordinates, a) is the reduced fre- 
quency, k the number of waves; L is a characteristic length, such 
as half maximum wing chord, for example. 

For harmonic oscillations 


0 = Lv exp (jvt + o)(3 2 x) £*(x,y,z) + Lvx (l8) 

p = p 0 v 2 exp (jvt + cop 2 x)p*(x,y,z) + p Q (19) 

£ = L exp (jvt + <cp2*)5*(x,y) + l£ 0 (x,y) (20) 

With these new variables, the Euler equations for harmonic oscil- 
lations become — 

(u> + J^jo*(x,y,z) + p*(x,y,z) = 0 


( 21 ) 



NACA OM 1402 


9 




( 22 ) 


Assume, for the present, an Infinitely thin lifting surface, lying 
approximately in the plane z = 0 , hence as having only infinitely small 
deflections £ . As boundary condition for the solution of equations (21) 
and ( 22 ), there is the given downvash on the lifting surface 


w 


dt 



(23) 


By equations ( 17 ), ( 18 ), ( 20 ), 
the lifting surface is 


and ( 23 ) the reduced downwash on 


w*(x,y) = 4*(x,y,0) 

Li 


7 ' 


1 - P‘ 


:|c> + |^ 5 * ( x,y) 


(2lf) 


The solutions of equations (21) and (22) are usually represented 
in integral form by means of Green’s functions or difference kernels K 

( 25 ) 

(2 6) 

Hie integrals are extended over the lifting surfaces; w*(x',y', 0 ) 
is the given reduced downwash on the lifting surface; p*(x' ,y 1 ,+©) is 
the reduced pressure on the upper side of the lifting surface. The 
pressure p*(x' ,y ' ,-0) on the bottom side is inversely equal. Hie 
Kernels and Kq must satisfy the wave equation (22), as is readily 

apparent. Let F(x,y,z) be a function that satisfies the wave equation 


®*(x,y,z) = / / K-l ( x - x',y - y , ,z)p*(x’,y , ,+ 0 )dx’ dy* 


P*(x,y,z) = // K 2 (x - x’,y - y’,z)w*(x , ,y',0)dx’ dy' 


(a 2 + K 2 )F(x,y,z) = 0 


(27) 
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and disappears at infinity. Ely equation (17) the reduced coordinates 
are 


y = -jy z = -jz real, if p > 1 


(28) 


Then the Kernel K^ is generally given by 

r x * 

K^x^z) = / exp a>(a - x)da ^ F(a,y,z) 
^ -00 

. . exp (-Jk \/x 2 + y 2 + z4 

F(x,y,z) = S= - L , if P«“l 

2rt \/x 2 + y 2 + z 2 

( / 2 "3-”?) 

F(x ’ y ’ 2 > ; f _a g ’ “ * > 1 

« \/x - y - z 


(29) 


(30) 


(31) 


Integration of equations (50 ) and (51) , with respect to y and y 
from -00 to +00 gives the corresponding functions for plane flow 


F(x,z) = - £ H 0 (2) ^ x/x 2 + z 2 ), if p < 1 
F(x,z) = -JJ 0 (k\/x 2 - z 2 ), if p > 1 


(32) 

(33) 


where Jq, Eq are cylindrical functions of zero order. 
The surface equation 


o-2-2 

x^ - y - z =0 


(34) 


gives a cone envelope, the so-called Mach cone. As the. roots in equa- 
tions (51) and (53) must always remain real, the range of integration 
of equation (25) in the supersonic range p > 1 must be limited to the 
Mach cone. The Kernel K>> in the supersonic range is 
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when 0 > 1 


n/x- - y 

By equation (35) > integration with, respect to y in plane flow 


(35) 


Kp(x,z) = -J jcn + Jq (k \fx? - z 2 ) when p > 1 


(36) 


3h the subsonic range p < 1, no simple formula in Cartesian coor- 
dinates can be given for Kernel Kp, since it is impossible to satisfy 
the Kutta wake condition in simple manner by this procedure. Curvi- 
linear orthogonal coordinates (u,v,w) must be introduced, making the 
lifting surface itself a member of the family of orthogonal surfaces. 

Let 


x = x(u,v,w) y = y(u,v,w) z = z(u,v,w) 
The auxiliary functions 

Xj-2 - fex\ 2 


5u 


(^A 2 , (bz 

\8u/ ydu 


V-2 = M 2 + M 2 + (%z \ 2 

\bv) ySvy 

2 = (bx \ 2 (by \ 2 l ^\ 2 

\d\rj \5uy \dv/ 


sire defined. 


(37) 


(38) 


For the three-axial ellipsoid with the half axes a, b, c we 
get, for example 


[P(u) - ej [ E (v) - ej E (w) - ej 


y 2 = 


( e l ~ e 2) ( e l " e 3) 

E(u) - e^J [g_(v) - e^j |] E (w) - e^ 

( e 2 

[g.( u ) - e^J 

" e l) i e 2 

|]p(v) - e^ 

e 3) 

[p(w) - e 3 ] 


( e 3 “ e l)( e 3 " e 2 j 


(39) 
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3®1 " 

b 2 + c 

/- 

CM 

1 

CM 



5e 2 = 

C 2 + 8 

, 2 - 2b 2 

► 


5e 5 - 

a 2 + b 

, 2 - 2C 2 


IT 2 = 

[p(u) 

- p(vjj 

£p(u) - 

p(")J 

v- 2 = 

[p(v) 

- p(v) 

'p(v) - 

g(uj 

w 2 = 

[p(v) 

- p(u) 

p(w) - 

p(vj 




(4l) 


£ is the elliptic function of Weierstrass. For the elliptic lifting 
surface c = 0. In orthogonal coordinates, wave equation (27) becomes 



In the cases involved here, this equation can be resolved by sepa- 
ration of the variables 


F = R^u^fvjw) 


(45) 


The coordinate u = constant is to denote a convex surface, u = 0 
the lifting surface, upper and lower side, R n (u) a retarded solution 
of the wave equation (42), which disappears for u = 00 . The functions 
S n (v,w) are then so-called surface harmonics of the surface u = constant j 

they satisfy the orthogonality relation 



dv dw 


f 0, if n ^ m 
k n , if n = m 


(44) 


The plane element is 
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Hie Integral can be extended over an arbitrary surface u = constant, 
because U/VW is independent of u, if tbe separation of the variables 
is possible. (Cf. eg.. (4l) . ) The surface coordinates v,w shall be 
unequivocal at least in certain ranges, as for example, the geographic 
coordinates of the sphere in the ranges 


0 < v < 2rt, - ^ < w < + ^ 


The coordinate w = 0 denotes the plane of symmetry, normal to the 
lifting surface, parallel to the flight direction, the equatorial plane, 
the coordinate u = v = 0 the leading edge and u = 0, v = it the 
trailing edge of the lifting surface. The elliptic coordinates given 
in equation (39) still are ambiguous, but can be made unequivocal by 
introduction of the Weierstrass elliptic c functions. 

Introducing, for abbreviation, the differential operators 


D x = 

% = 


5 5 ^ 

dx dx* 

,5 

JL +. _d_ 
By By' J 


(^5) 


Then there follows for the difference kernels and K£ of the 

integrals ( 25 ), ( 2 6 ) the differential equations 


DxK = DyK = 0 


(W) 


To satisfy these conditions in curvilinear orthogonal coordinates 
too, the characteristic function 


RqCu) 

G(u,v,wjv',w*) = 2_ rT(o} SnOv^Sj^v' 




(47) 


is introduced. 

By definition this function satisfies wave equation (42)} it further 
shall satisfy the compatibility condition 
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C^G — CgDyG — 


G v .(0,Jt,0jv’,w* ), 
M 0,0,Ojv',w'), 


Gy.(u,v,wj 0,0) 
G v » (u,v,wjit,0) 


(48) 


Ey appropriate normalization of G the constant can be made 

equal to unity . Then the Kernel Kg in the subsonic range 3 < 1 can 
be represented as follows. (Cf. ref. 16.) 


K S *■ - gjrjQCUjTjVJTSw') + 

<V (u,v,wj0,0) |G v .( 0 ,jc,0;v , ,w , )T 1 (K,a)) - Gv-(0,0,0jv' ,v* )] 
Kg = -Agj^o - ^ 7 ^G(u,v,w;v , ,w') + 

Gy, (u,v,wjrt,0) [(h(0,jt,0!v',w') - G v .(0,0,0jv',w , )T 2 (k:,cd)| 


> (49) 


o-jjt/2 


Tl ( K ^) = 


f _ exp ok(u, 0,0)0^. (u,0,0;rt,0)^^^- du 


U(u,0,0) 


r“- ifl/2 -v / u q o) 

/ exp uk ( u,0,0)G^ nr >(u, 0,0; 0,0)- - *- ■ ' ■ 

^-ooHat/2 U(u,0,0) 


du 


(5P) 


T 2 (k,cu) = 


-v<»-ijt/2 V(u,0,0) 

exp ok ( u^O^ 0)G^^ i (u^rt ^ 0; 0; 0)--p tt &U 
-oofijc/2 U(u,0,0) 

' oo jLyc /2 ” 

exp uk(u > 0 > 0)G vv i (u,3t,0}Jt,0) ^ - ^^| du 
-«*irt/2 • U(u,0,0) 


(51) 


These equations hold for complex values of v, k , a) corresponding 
to equation (l6), i . e . , for damped and amplified oscillations also. 

For lifting surfaces having a symmetry axis x = z = 0 , for the 
elliptic lifting surface, for example, Tj_ = -Tg. For the constants A-^ 

and Ag in equation (49), the following condition is applicable 
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- — 7—— = // dx dy lim G(u,v,w;v' ,w' ) 
^•1 + ^2 z=0 ° z 


= JJ Gyx(°> r > v > yX >*') ^ dv dv 

The Integrand must he a Dirac 6 function 
G u (0,v,w;v’,w') = 


(52) 


if 

V ^ V* 

or w ^ w * 

if 

V = v' 

and w = w 


(55) 


The conditions (52), (55) are such that the Kernel closely related 
with K 2 of the integral representation of the velocity potential as 
function of the downwash always reproduces the given downwash on the 
lifting surface . (Cf . ref . 15 • ) If both the leading and trailing edge 
of the lifting surface are in uniform flow, we get 


A x = Ag 


(5*0 


hence zero lift in steady flow. Unsteady motion of lifting surface 
creates additional members, among them the aforementioned Kelvin impulses, 
which are acceleration forces of the comoving fluid. The solution (5*0 
is termed Kelvin solution; in the 19th century it was regarded as a para- 
dox of hydrodynamics . Experiments prove that the Kelvin flow occurs at 
the start of the motion, but is soon changed by the boundary layer of 
the lifting surface. Kutta therefore made the phenomenological assumption 
of smooth flowoff at the trailing edge, which leads to the constant 

A 2 - 0 ( 55 ) 

Kutta' s solution (55) gives a 20- to 50-percent higher circulation 
of the lifting surface; every steady and unsteady lifting surface theory 
up to now rests on this solution. However, actually the trailing edge 
is in a weak flow, so that 

A 1 > a 2 > 0 ( 56 ) 

A second conditional equation for the constants A^_ and A 2 of 
the general solution (*4-9) might be gained either from an unsteady 
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"boundary-layer theory or from measurements of the unsteady pressure dis- 
tribution. But the viscosity of the fluid discounted for the present 
makes itself felt considerably in the final result. 


5. TWO -DIMENSIONAL SOLUTIONS 


The wave functions Rq and Sn are not veil known except for the 
case of the strip and the circular plate. 

For the strip, there are the coordinates of the elliptic cylinder 


x = ‘■cosh u cos v y = w z = sinh u sin v] 
U -2 = V 2 — sinb^u + sin^, ft = 1 


( 57 ) 


With these coordinates, the characteristic function for two- 
dimensional flow reads 


00 Ue n (2 )( r ,u) 

G(u,v,v') = 2_ ToTTT"* — C Be n( K > v ) se n( K J v ' )> if K > 0 

n=l Nent ' (k,0 ) 

o(u,t,t-) - i m 22% u - °° s < T - v ]\ , K ,. 0 

7 7 4 cosh u - cos (v + v 1 ) 9 


( 58 ) 

( 59 ) 


Ne n and Se n are Mathieu functions in the Goldstein norm. McLachlan 
employs the parameter q. = k 2 /4. The tables of the NBS (ref. 17 ) contain 
these functions in a different norm with the parameter s = k 2 . 

By equations (2 6), (50), and (58), the Kutta condition (55) gives the 
reduced pressure. (Cf. ref. 14.) 


p*(u,v) = -- cf w*(0,v' )sin v' dv' { (00 -. — ^-)g(u,v,v') + 
rt Jo I \ sin v’ dv'/ 

( u > v >0) (Gv(0,jt,v' )T(k , 00 ) - Gy(0,0,v ' )j 1 


(60) 
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IT 


T(k,03) 


poo-irt/2 


pOO_i]t/2 

-aH-irt/2 


exp (-03 cosh. ^G^i (u,0,it)du 


exp (-to cosh u)G vv t (u,0,0)du 


(61) 


For k = 0, i.e., for incompressible fluid, the characteristic 
function (59) is inserted in equation ( 60 ) , so that the solution becomes, 
term by term, the Kussner- Schwarz solution. (Cf . refs. 1 and 14. ) For 
k = 0, there are large five-place tables of the aerodynamic coefficients 
of the wing with aileron and tab on the assumption of the bent flat plate 
as substitute system. (Cf. ref. 11.) For k > 0, Blanch and Fettis 
(ref. 18), as well as Timman, van de Vooren, and Greidanus (ref. 19), have 
published five-place tables of coefficients for the flat plate and 
NLL Amsterdam (ref. 19a) also for the bent flat plate. 

For small aspect-ratio wings, such as tapered delta wings, the speed 
and pressure changes in x direction can be approximately disregarded 
and so eliminate compliance with the flowoff condition. By utilizing 
the orthogonal coordinates 


y = cosh u cos v z = sinh u sin v ( 62 ) 

the Kernel becomes then approximately equal to the regular solution 

of the wave equation 


K 2 (u,v,v') = -— G(u,v,v r ) ( 65 ) 

where G is the function given in equation ( 58 ). Equation ( 65 ) holds 
for all Mach numbers p. The lifting surface with its given downwash 
distribution is divided in separate strips parallel to the y axis; 
each strip being treated according to equations ( 26 ) and ( 65 ). Corre- 
sponding solutions and tables have been computed by Merbt anfl TjmriAhi 
(ref. 20). How far the approximation ( 63 ), which equation (46) does no 
longer satisfy, is practicable must be left to exact solutions. 


6. THREE -DIMENS IONAL SOLUTIONS 


For k = 0, the wave equation reduces to Laplace's potential equa- 
tion. For the circular plate, the solutions are spherical functions and 
for the elliptic plate the Lame functions, so that corresponding charac- 
teristic functions can be set up. For the circular plate, Schade (ref. 21) 
computed solutions with spherical functions by a complicated procedure. 
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When k > 0, solutions for the circular plate can be computed with 
spheroid functions developed by Meixner. For the elliptic plate, which 
is of greater aeronautical interest, the corresponding wave functions 
are still lacking. 

Kussner (ref. 21a), Reissner (refs. 22 and 23), and Dengler -Goland 
(ref. 23a) have published approximate solutions patterned after Prandtl ' s 
vortex filament theory for high aspect-ratio wings. The closed solution 
of the two-dimensional problem of airfoil theory is assumed known and a 
correction of this solution is computed by resolving one or more linear 
integral equations. When the Wronsky determinant of the cylindrical 
functions is taken into consideration, the theories of Kussner (ref. 21a) 
and Reissner (ref. 22) are identical in every respect. For incompressible 
fluid k = 0, the method is tested. But for k > 0 it is very complicated 
and has not been carried through numerically up to now. Compressibility 
effects, in aeroelastic problems, do not make themselves felt appreciably 
until P > 0.7 (ref. 3). This is the reason for the practical results 
obtained thus far with the wing-flutter theory at the pressures in the 
range 0 < {3 < 0.7 computed for k = 0. 

In the supersonic range p > 1, the Kernel K2 is known in the simple 
form (eqs. (35) and (36) ); hence the many tables of the aerodynamic coef- 
ficients for plane flow available now (cf. Weber, refs. 23b and 5)> as well 
as three-dimensional solutions for rectangular, triangular, and swept - 
wing conf igurations . (Cf . Watkins (ref. 8); Lomax, Heaslet, and Fuller 
(ref. 24); Kelson (ref. 25); Watkins and Berman (ref. 26); Miles (ref. 27); 
Walsh, Zartarian, and Voss (ref. 27a).) If the flutter study involves 
thin plates - a problem particularly posed in the supersonic range - the 
integrations of equations (26) and (35) must be made from case to case, 
corresponding to the wing contour and the number of points of the elasto- 
mechanical deformation matrix. This is apparent from the flutter theory 
developed further on. Solutions for P = 1 are most easily obtained from 
the supersonic solution by the boundary transition p — >1. (Cf. refs. 28, 

29, and 30. ) 


7. IMPROVEMENT AKD EXPERIMENTAL CHECK OF THEORY 


In incompressible fluid (k = 0) and two-dimensional flow, profiles 
of finite thickness can be treated as unsteady by the conventional methods 
of mapping (refs. 31 and 32). Numerical solutions for small oscillations 
have been computed by Couchet (ref. 33). The hope of obtaining aerody- 
namic forces that are in better agreement with experiments has not been 
fulfilled at all, or only in a small part. (Cf. ref. 34.) In fact, if 
a frictionless airfoil theory is to be maintained at all, it seems 
advisable to relinquish the Kutta condition. A possible improvement of 
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theory by the general solution (^ 9 ) with the phenomenological con- 
stant Ag/A-^ has been pointed out previously. 

Drescher (ref. 35 ) measured the pressure distribution of an oscil- 
lating wing with control surface in plane flow in the water tunnel and 
compared it with lifting surface theory at k = 0. The amount and the 
phase angle of the differences vary very little in the reduced frequencies 
range oi = 0 to 2 . 5 j this means that the differences correspond to those 
in steady flow. Similar results for total lift and moment of a rigid 
oscillating airfoil have been obtained by Greidanus, van de Vooren, and 
Bergh (ref. 36) in wind-tunnel tests. The Kelvin impulses were measured 
at v = 0 ; they are in fairly exact agreement with theory, up to a small 
additional damping portion. Dtirr (refs. 3 ^ and 37 ) studied the flutter 
of a wing with and without aileron in plane flow in the wind tunnel and 
found' agreement in several cases between computed and measured critical 
speeds after multiplying the theoretical pressures by the reduction 
factor = 0 . 7 . However, this agreement is contingent upon the choice 
of mean camber line used as theoretical basis for the aileron with 
internal ba l ance . The experimental findings may be attributable to 
profile shape effects or to boundary- layer effects. A new theory which 
accounts for both is therefore desirable. 


8. MATRIX THEORY OF WING FLUTTER 


The earlier flutter theory haB been described in several compre- 
hensive reports (refs. 1 , 4 and 5 ); they chiefly rest on the replacement 
of the elasto-mechanical system of the wing by an elastic beam and on 
the assumption of plane flow past the individual wing sections. It 
results in a system of linear differential equations that can be solved 
with any desired degree of accuracy by iteration. These simplifying 
assumptions suggested themselves as a result of idle wing design of the 
time on the one hand and the absence of large analog computers on the 
other. The demands of flutter theory are one of the driving forces in 
the evolution of the first automatic computers in Germany and the U.S.A., 
of which a considerable number is now available. This circumstance has 
led to a marked improvement of the elasto-mechanical substitution system 
a n d of the theory of natural oscillations of the airplane. And here is 
where the matrix calculation is particularly suited. The first step 
that must be taken is to improve the aerodynamic principles of the flutter 
theory. The aerodynamic forces which must be inserted in the principal 
equation ( 7 ) in order to be able to resolve aeroelastic problems, must 
also be represented in matrix form. This is theoretically possible, as 
soon as the Kernel K 2 in equation ( 26 ) is known. 
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To this end, the deflection of the lifting surface must be repre- 
sented in appropriate manner. Assume a given class of functions 


F(x,y; r,s) 


(64) 


vhich are regular in the range of -the lifting surface, i.e., finite, 
continuous and differentiable; x and y are the reduced coordinates 
introduced in equation ( 17 ) ; r and g are integral parameters, starting 
with r = s = 0. The reduced displacement £* of the lifting surface 
is developed with respect to a series of functions F by means of the 
formula 


5 *(*,y) = F ( x >yJ r i> s i)Sik££? + a ± < m (65) 

are the given displacements in n 

of the lifting surface. The functions 
so chosen that their determinant 

det|F(x k ,y k jri,s. L )| ^0 ( 65 a) 

By equation ( 65 ), the insertion of the coordinates x = x^, y = y^ 
gives the relation 


= ( m 2 1 ) ciLOSe:n points ^ x k ,y k ^ 

F and the point coordinates are 


F ( x ZW r i> s i)gik = 8 Zk 

and with the abbreviations set up in equations (3) and (4) the matrix 

g in * [ p (w r i' e 41 - 1 < S6 > 

The representations ( 65 ) and (66) are exact, when the function £*(x,y) 
belongs to class (64) and of lower than mth order; otherwise, it is 
approximately valid. When the lifting surface displacement £* contains 
a discontinuity due to a bend in the control surface, two separate rela- 
tions must be established and then later combined along the control 
surface axis. 
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Integration of equation ( 65 ) with, respect to the lifting surface 
gives the well hnovn numerical integration formula 


dy = ^*(x k ,y k ) = ( 67 ) 


with the weight factors 


a k = s ik JJ F ( x >yJ r i> s i)dx <3y (68) 

When f (x,y) = £*(x,y) is an arbitrary function of the class (64) 
and of lower than mth order, equation ( 67 ) is exactly valid, otherwise 
approximately, the error can be computed by Mises ’ method (ref. 38). 

In numerical integrations polynomials, i.e., linear forms of the 
function 


F(x,yjr,s) = xFy® 0 ^ r + s < m ( 69 ) 


are generally employed . When the lifting surface u = 0 is given in 
curvilinear orthogonal coordinates, it is, however, more advantageous to 
utilize linear forms of the function class 


F(x,yjr k ,s k ) = U(0,v,w)%(v,w) 


(TO) 


in which are the previously introduced surface harmonics. These 

functions too are regular in the entire lifting surface range. 

Assume an arbitrary integrable mass distribution m(x,y) over the 
lifting surface. 5his is to be replaced by n single masses Mjj. in 
the points ( x k ^y k ) i* 1 such a way that the total masses, static moments, 

moments of inertia, etc. of both distributions correspond. This condi- 
tion is met when 


MfcF ( x k>ykJ r i> s i) = j(/" F ( x ,y* r i ,s 1 )m(x,y)dx dy (71) 
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M k 


■ Si k JJ'f (x,y}v lt B£m.(x,y)Qx dy 


(72) 


Equation ( 72 ) is similar to equation ( 68 ) . 3h contrast to the 
aforementioned function f (x,y) which must he regular,, function m(x,y) 
in equations ( 71 ) and ( 72 ) may also contain integrable singularities 
(single masses). These integrals are called Stieltjes^ integrals. When 
the reduced coordinates equation ( 17 ) are used, the right-hand side of 

equation (72) must he multiplied hy L 2 [\/l - p 2 , in order to obtain 
physical masses on the left. 

The same procedure is used for the replacement of the pressure 
distribution by n single forces. Assuming an aerodynamic replacement 
force of 


P k = exp fovt + o£ 2 x( k )jp* 


(73) 


equations (19) , ( 72 ), and (73) give the reduced replacement force 

2o v 2 ! 2 PP 

P* = g g lk jfjE^x,y;r 1 ,s 1 ^p*(x,y, + 0)dx dy (74) 

Insertion of equations (24), ( 26 ), and ( 65 ) in equation (74) gives 
p n v 2 L 2 p n v 2 L 2 

- f— 5 %lt (T5) 

X - p X - p 

h±Z = 2 Jf W [ X, 7 }V i ,e i) &X ^ *Jf (“ + i x ' ’ 7 ' dy ' X 

Kg(x - x',y - y',0) ( 76 ) 


The factor 2 in equations ( 74 ) and ( 76 ) is based on the fact that 
the integration covers only the upper side (z = +0) of the lifting Biirface, 
while the pressures on both sides of the lifting surface are inversely 
equal. In the supersonic range 3 > 1, the Kernel Kg is always known 

and given by equations ( 35 ) and ( 36 ) in cartesian coordinates x and y. 
Therefore, it is best to use polynomials in x and y as development 
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functions. (Cf. eq. ( 69 ).) In. the subsonic range p < 1 , the problem, 
is more difficult. Curvilinear orthogonal coordinates (37) function 
class ( 70 ) must be used as development functions. Partial integration 
of equation ( 7 6 ) gives the aerodynamic matrix 

\i =J] r s i (v,w)^dvdvx JJ S l (v*,w«)^av« dw« X 

^5 - ^r^Co^wjvSw*) (77) 

Due to the insertion of the Kernel Kg according to equation (49) in 
equation ( 77 ), all integrations can be carried out in closed form by- 
means of the orthogonality relations (44). The integrations are extended 
over the entire orthogonal surface u = 0 , i.e. , both sides of the lifting 
surface . Thus, matrix h^ can be computed from a limited number of 1% 

values, according to equation (44), and of the wave function values 
appearing in equation (49) . In the case of a hinged aileron, hence unsteady 
deflection £*, the integrations over both areas of the lifting surface 
must be made separately. Equation (44) is then no longer applicable and 
infinite series are obtained. Only in plane flow can those integrations 
be carried out in closed form. 

If the Kernel Kg is not known, the integrodifferential equation 
v*(x,y) = y.h--.-. L 4- J^k*(x,y) 

V 1 _ p2 \ x / 

- llm “ ff K l( x - x '>y “ y*,z) P*(x,y> + 0 )dx‘ dy‘ ( 78 ) 
z =0 oz 10 


following from equations (24) and ( 25 ), must be resolved, where K^ is 
given by equations ( 29 ) and (30) 

Integrating while using equation ( 27 ) gives 


£*(x,y) 


= i/ K o (x ~ 


x’,y - y')p*(x , ,y‘, + 0)dx* dy* 


(79) 
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%(x,y) 




da exp co(a 



5 2 _ a 2 ^ exp (-J* y4 2 + y 2 ) 

aa 2 ay 2 / ^2 ~ y 2 


( 80 ) 


If the Kernel were regular everywhere, it could he used as 

development function P and equations (74), (75) > and (79) would then 
immediately yield the aerodynamic matrix 



- h) 


-1 


But Kernel Ko(x,y) has a singularity at x = y * 0. The reverse 
problem must therefore be solved some other way, which cannot be 
discussed here. 

The sum of the reduced mass and aerodynamic forces is, according 
to equations (6), (25), (75), and (75) 


P* = -IM 
k 


(k) v2? £ 




P Q v 2 L 2 
1 - P 2 




(®L) 


On the other hand, elas to -mechanically, we get by equations (2), 
(25), and (73) 


«f - b ik pk 

(82) 

( ik) exp - xj 

(83) 


Multiplying equation (8l) by b_j^Jv 2 L and inserting it in equa- 
tion (82) gives the principal equation of the flutter theory 




<») 



MCA HM 1402 


25 


iF 


Pq v2l 

v 2 (l - p 2 ) 



(l - P 2 ) 


(85) 


The parameter p has the dimension of a mass. When p = 0, equa- 
tion (84-) becomes equation ( 7 ) • Equation (84-) is the desired matrix 
equation for defining the n eigen values 7^, of the oscillating wing. 


9. ANALYTICAL SOLUTION METHOD 


The bracketed matrices in equation (84-) can be computed for given 
values of k amj co and are complex. Hence, it seems logical to ’ 

consider equation (84-) as a complex eigen-value problem, whereby the 
aerodynamic matrix is not self -adjoined. Prior to 194-5 , little was known 
about complex eigen-value problems. Since then, Wielandt (ref. 39) has 
developed appropriate iteration methods for calculating complex eigen 
values and proved their convergence. Recently, Gossard (ref. 4-0) applied 
this procedure and proved it again. 

Equation (84-) thus yields a set of complex eigen values whose 
azimuth angle is the assumedly known damping phase angle gj_ 


exp jgj^ 



( 86 ) 


These physical solutions are obtained from the group of mathematical 
solutions by graphical interpolation, 7 ^(k, 0)) being plotted in the 

complex numerical plane. The critical speed v-j_ is computed from the 
obtained values v^, tc, and ca, according to equation (l6). The check 
on whether 


P 


_ jK 
Cq 0) 


( 87 ) 


conforms with the assumption is made by graphical, method. And, since 
the airplane is to fly at different altitudes, the air density Pq must 

be varied too. 

The conditions (86) and ( 87 ) confine the physically possible solutions 
materially, so that often no physical solution is found in the explored 
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speed range. Ibis is exactly -what the aeronautical engineer wants. To 
be on the safe side, however, he would still have to vary the most impor- 
tant parameters of the airplane with their limits of error and manufac- 
turing tolerances. This can.be done by one of the customary disturbance 
calculations, without having to repeat the whole calculation (ref. 4). 

It might be noted that Wielandt (ref. 4l) has developed mathematical 
methods for estimating the upper and lower limits of the complex eigen 
values 7vj_, when the principal equation (84) is given. When these limits 
meet the technical requirements, the solution of the characteristic value 
problem can be omitted. 

Bearing in mind the structural damping, the reduced elastic force 
follows from equations (l), (20), and (73) at 

_• 

P* = I£fa(ik) expja3p 2 (xi - x k ) + jg] = L^a.^ (88) 


The sum of the elastic, mass and aerodynamic forces is then, 
according to equations (8l) and (88) 

2 Po v2l2 

- “oo v - rri? ' “5 (89 > 


Flutter is generally initiated by small external forces. Hence, 
the thought lies close to introduce a given periodic outside force Q* 

!*, for example, 

or a force distribution representing a given "sinusoidal" gust. The 
amplitudes £*(v,v) of the lifting surfaces can then be computed as 

functions of the frequency and speed from the inhomogeneous equation 
system ( 89 ). 

Near the critical flight attitude, the amplitudes increase consid- 
erably and become infinite in the flutter range by reason of the disap- 
pearance of the denominator determinant of equation ( 89 ). This method 
also affords a complete survey on the flutter behavior of an airplane. 
This method is especially favored in Russia. It may also be used for 
investigating the effect of atmospheric turbulence, by analyzing it 
harmonically and posting it on the right-hand side of equation ( 89 ) . 

The harmonic components of the elastic stresses calculated from the wing 
deflections are statistically superposed. 

The six degrees of freedom of the aircraft motion as a rigid body as 
well as the free control surface notations must also be borne in mind in 


as interference function in equation (89), such as Q* 
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every flutter study. Up to now, the theory of path motions of airplanes 
rested on these degrees of freedom. However, no sharp boundaries between 
flight stability, flutter, and gust stresses can be drawn; there are many 
transition modes. Collar (ref. 42) said: "in short, it is evident that 

we are no longer dealing with a series of subjects each in its own water- 
tight compartment: there is a definite coalescence of the subjects into 

an integrated whole, which may be defined as the dynamics of a deformable 
airplane. And we are faced with the question: what are to be our methods 
of treatment of this unified problem? " 

Rea (ref. 43) demonstrated how by Laplacian transform all three 
do m ains of aeroelasticity can be jointly treated (Transfer Function- 
Fourier Method) . The transfer function of a dynamic system transforms 
the input into output, for example, control hinge moment into tail sur- 
face force. For the numerical solution of these and, in fact, of all 
linearized problems, the use of electrical analogies is recommended. 

They rest on the formal similarity of the differential equations of an 
ordinary linear mechanical system with a finite number of degrees of 
freedom (Lagrange’s equations) and the differential equations of a linear 
electric network (Kirchhoff’s law). The displacements represent the 
electric voltages, the generalized forces, and the electric currents. 

For conservative systems, the analogous electrical network contains only 
passive elements, resistance, capacitance, inductance, and transformer. 
When air forces are involved too, amplifiers must be included. 

The solution, in characteristic value problems of flutter theory, 
is obtained by connecting periodic interference currents Qg with the 

network and measuring the currents and voltages. One particular advan- 
tage of such analogy devices is that structural changes in the network 
can be easily copied in the design stage of the aircraft. The error of 
measurement increases nearly with the root of the degrees of freedom, 
therefore, is still conservative even for 100 degrees of freedom. McNeal, 
McCann, and Wilts (ref. 44) described the analogy setup developed at the 
California Institute of Technology -which carried out flutter calculations 
with 70 degrees of freedom. 

The stability of the solutions of linear differential equations is 
defined by Routh’s criterion; but it becomes inconvenient with a large 
number of degrees of freedom. Nyquist, therefore, subjected the differ- 
ential equations to a Laplacian transform and developed its characteristic 
equation in powers of Laplacian operators. The stability of the solutions 
of the particular system can then be proved by a simple theoretical func- 
tion operation. Nyquist ’b method has been used for some time in electro- 
technics, and has been applied also to flutter calculations by Dugundji 
(ref. 45). 
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10. FLUTTER OF SKIN PANELS (COVERING) 


The covering of an airplane usually consists of thin sheet panels 
attached to frames; they can perform aeroelastic oscillations even on 
a rigid frame. To the extent that the problem can be linearized and the 
wing covering is approximately flat, this type of flutter could also be 
treated by the matrix method developed in section 8. However, this would 
require the assumption of a rather large number of degrees of freedom 
(points), and so the solution of very great matrices. 


For this reason we examine a simpler, ideal case, namely, the infi- 
nitely flat thin elastic plate, pin supported on an infinite rectangular 
lattice framework with the spacings and 1 2 - The mean normal 

stresses Tp/h and T^h and the plate shearing stress S/h may be 
regarded as constant for very small deflections £, even if the lattice 
framework is immovable. The plate is of h thickness, and M is the 
mass per unit area. The linear differential equation Jref . 46) for £ 
reads then 


a 2 ? 


4 d 2 f; 

B7 £ + Ti — 7 - + 2S 

8x 2 8x8y 


+ T 2 — §■ + M — §• = -H(x,y,t) 
2 8y 2 St 2 


(90) 


x, y, and z are ordinary cartesian coordinates; z and £ are counted 
positive upwards. Compressive stresses are computed positive. The 
bending stiffness 


mG h? 

B ’ T Ig 


(») 


is assumed complex for periodic processes, to allow for structural 
damping. The pressure jump on the plate is 


00 r, 

H(x,y,t) = p(x,y, + 0,t) - p(x,y, - 0,t) + J 6 (x - nlj <^(y,t) + 

n=~"oo v ; 


r 5 (y " mZ 2)Sn^ x ^ t ) 


(92) 
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p = aerodynamic pressure, 5(x) Dirac's 5 function, (^ n (x,t) the 
normal force per unit length exerted, on the mth "bullhead, (longitudinal 
member) of the lattice framework. In the subsequently involved periodic 
solutions, the functions and repeat after each two steps. 

Equation (90) is solved by the substitution 


00 oo 

(x,y,t) = ReRe’ ^ A rg exp (ivt + irkjX + jak^ 

—00 —00 ' / 


i 2 - ±i = i 2 = -i 


(93) 


Apg are hyper complex constants; in addition 


\ = «Al ^2 = «/ l 2 


Continual support of the plate on the lattice frame is contingent 


upon 


00 CO CO 00 

—00 —00 / _on m 


— 00 —CO 


A 2r+l,2s 0 


00 00 00 00 

'L- ^ A 2r+l,2s+l " ° 




The ind ividual members of (93) represent traveling waves moving in 
x direction at speed -v/rk^. The air strikes the upper side of the 
plate at constant speed v along the positive x axis, while still air 
prevails at the lover side. A new system of coordinates 


= x + 


vt 

rk]_ 


x 


y = y 


z = z 


(95) 
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that moves along with the traveling wave r is introduced. In the over- 
scored coordinate system, a steady flow problem exists, the upper side of 
the static wave r is in streamwise flow at speed 


v + 



the lower side at speed 



For k = 0, equation (22) gives the differential equation of the velocity 
’ potential 





4>(x,y,z) 


= 0 


(96) 


from equation (25) follows the boundary condition for z = 0 


w(x,y) = $z(x,y,0) = %y?(x,y) (97) 

By equations (ll) and (15), the pressure on the upper _ s id e of the plate is 

p(x,y, + 0) = -p 0 v$ ? (x,y, + 0) + p Q (98) 


The particular deflection of the plate is assumed as 


£ (x,y) = exp (ir^x + jsk 2 y) 


(99) 


and the corresponding velocity potential for z 0 as 


. iEL. s 

®(x,y,z) = ■ - U. -7 — exp 
p 0 vrlci 


irk^x + jsk^y - z //l - ^\r 2 k^ + s 2 k| 


( 100 ) 
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As a result of equations ( 96 ) to (100) , tlae particular pressure is then 


p(x,y, + 0) = EpgCviexp^k-LX + 



(101) 


with 


B (v) 
rs v 1 



+ s\| 


(102) 


when the ra&icand is X), and 



(103) 


when the radieand is <0. 

Putting equations (92), (93), (95), (99) , and (101) in equation (90), 
multiplying by 

exp^-irk-jX - Jsk^y) 

and integrating with respect to x,y from -» to +“, yields the 
coefficients 

^ - (Oo + (-DX + (-D 6 0 2 ] (Hj" 1 (104) 

■^rs = (A| ♦ .*4)h - A&. - 2rsk-i_k2S - 

‘S&S * + Brs (105) 


The hypercomplex constants determine the type of the solutions. 
Corresponding to the support conditions (9*0, there are four different 
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types of form changes, hence four types of flutter. Only the straight 
oscillation modes are considered here. By equations (9*0 and (104) 



( 106 ) 


Energy absorption from the. airstream and flutter is possible only 
when at least one of the coefficients B^g becomes imaginary. (Cf 

eq. (103);) Equation (106) then becomes complex, from which follow two real 
equations for the determination of the eigen values v and v. For 
v = g = 0, only static instability is possible in the subsonic range, 
but no flutter. By equation (102) all coefficients B^g are then real, 

hence the sole result is a real equation (106) for computing the critical 
speed v. "When at least two f met ions N rs of the same parity disappear 

simultaneously, equation (106) itself is satisfied, since » _ 00 can be 
any number. 

Hedgepeth, Budiansky, and Leonard (ref. 4-7) developed a representa- 
tive theory for plane flow (k 2 = 0) and computed the stability range of 
the fundamental oscillations r = -2,0,2 for a large number of parametric 
values on the basis of the Nyquist diagram (ref. 45). !Ehey found that 
structural damping g could lower the stability range under certain 
conditions. The limitation to fundamental oscillations for reasons of 
simplicity is a first approximation. [Hie support reactions of the lattice 
frame cause discontinuities of "the higher derivatives of the deformation 
area £(x,y,t), to whose representation Fourier series with infinitely 
many terms corresponding to (93) are required theoretically. Since N rs 

is a 4th order polynomial in r and s, the series (106) converge fairly 
rapidly. 


When external air forces T > 0 act on the plate or it is heated 
by skin friction, the plate can buckle statically even without air forces. 
By equation (90) the differential equation of the static plate (ref. 46) 
reads 


+ Ti 



2S 




Sy 2 


= 0 


(107) 
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Inserting equation (93) in ( 107 ) gives 


rs 


jf 


r 2 k^ + 


s 2 k| 




«2-i,2| 


; k£T x - ^sk^S - s 2 k§T 2 


= 0 


( 108 ) 


Finite amplitudes A rs are possible only vhen the bracketed term 
in equation (108) disappears, i.e., vhen T-j_ and T 2 have a critical 

magnitude. If these increas.e gradually, buckling vill follow at the 
smallest possible values of r and s, and a corresponding static wave 
field will be the result. Visualizing this wave field in a flow along 
x at speed v, its steady pressure distribution can be computed by 
equation ( 101 ), where x = x, v = v, v = 0 . The aerodynamic forces 
continue to deform the given wavy plate until a new state of equilibrium 
is reached or, failing that, until the waves change into opposite posi- 
tion (£ — ^ -£ ) • This loss of Btatic stability is a supersonic phenomenon 
which comes into being through the change in sign of the radicand in 
equation (103 ) . In the subsonic range the aerodynamic forces generally 
have a stabilizing effect on a static wave field, with exception of the 
aforementioned case of static instability, where a monotonic increase in 
wave height occurs until the linearized formula ( 90 ) becomes inapplicable. 

Fung (ref. i+8) investigated the static stability of a sinusoidal 
half wave in two-dimensional supersonic flow (k 2 = 0, v > Cq) ; however, 
his theory is identical with that of an infinite wave field. By assuming 
immovable supports and finite small amplitude, the problem becomes 
nonlinear with respect to the deflections, but still linear enough approx- 
imately as far as aerodynamic pressure is concerned. The critical speed 
of static instability is so much lower as the initial amplitude of the 
wave field is lower. Because the opposite position (t, — > -t, ) of the 
changed panel field has just as little stability, a complicated periodic 
process develops which is not harmonic and could not be calculated up 
to now. Such a process would lead to a rapid destruction of the shell 
(or skin) as observed on German V-2 rockets. Skin buckling on supersonic 
airplanes in all flight attitudes should therefore be prevented by appro- 
priate design measures. 

To the extent that the ideal case of the infinite plate with lattice 
frame investigated here is applicable to real aircraft, skin flutter can ' 
be considered as a typical supersonic phenomenon. Flutter processes can 
occur only below the critical speed of static instability. 
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11. STATIC OSCILLATION TEST 


The static oscillation test on airplanes is the oldest and most 
extensively employed experimental method of aeroelasticlty. Since 
flutter occurs mostly near static oscillation frequencies, the static 
oscillation test orig inally served to detect oscillation modes suspected of 
flutter, later on largely as experimental check of the design regulations I 
and II as a safeguard against flutter. Owing to the complicated structure 
of aircraft, there was no other possible way to meet these requirements. 

The measured static oscillation modes P^(x,y) were utilized occa- 
sionally as starting function for the flutter theory, i.e., the wing 
deflections t, were represented as linear modes of P^(x,y). The mass 

distribution m(x,y) of the wing is comparatively easy to define from 
the design drawing, from which the substitute masses in n selected 

points can be computed by equation (72). After inserting F^ 

Mj., and the measured natural frequencies in equations (l), ( 5 ), and 

(6), the stiffness matrix a^ can be computed without having to make 

stiffness measurements or calculations with their attendant sources of 
errors and approximations. To be on the safe side, every measured static 
oscillation mode from the first to the nth should be introduced as degrees 
of freedom, which appears possible by electronic computers. In order to 
be able to carry out this procedure, the measured static oscillation modes 
must, first, be corrected and orthogonalized (ref . 48a) . 

For the static oscillation tests, the readied airplane is suspended 
in a hanger from very soft springs or else mounted on a suitable elastic 
base and excited in one or more points by oscillators with slowly 
increasing frequency. At the natural frequencies Vj_, the ampli- 
tudes F^(P) are maximum. These can be measured or recorded successively 
and thus give a complete survey of the static oscillation modes F^(P) 
of the airplane, P denoting any one point of the airplane. 

The original oscillators were rotating unbalanced masses. But they 
were later replaced by electrodynamic oscillators, which are more accurate 
and easier to use. More recent oscillators are those developed by Herschel- 
Schweizerhof (ref. 49) and compared in an article by de Vries (ref. 50)* 
Originally, the amplitudes were read from so-called oscillometers mounted 
at several points of the airplane, which were replaced later on by 
commercial electrical pickup and recording instruments (ref. 5)* These 
pickups are either accelerometers operating with piezoelectric crystals, 
strain gages, or acceleration -responsive electron tubes; or else they are 
speedometers operating with induction coils. Accelerations and speeds 
are changed to deflections by electrical integration switches. Direct 
recording strain gages are suitable for elongation measurements. 
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These electrical Instruments speed up the natural oscillation mode 
measurements, increase their accuracy, and make phase angle measurements 
possible (ref. 51) • Only the natural oscillation excited in resonance 
has the phase position ~jt/2 + tot relative to excitation, while other 
accidentally excited osci lla tions have the phase position ~0 + tot and 
can be Interpreted as disturbances and eliminated. The exact measurement 
of the excitation energy also enables the damping phase angle g^ to be 

determined when the generalized mass of the oscillating system 

Pi = V F?(P)m(p) 

P 


is known. This can be obtained by application of known additional masses 
on the ’ oscillating system and by frequency measurements (refs. 52 and 52 a). 


12. MODEL TESTS 


When an experimental clarification of the most important aeroelastic 
problems is demanded before a new type of aircraft is f inished, a dynam- 
ically similar model of the planned airplane can be built. Its purpose 
may be the determination of the static oscillations, and can be met by 
the mathematical problem described in section 3. The outside of the 
model does not have to be geometrically similar to the full-scale design; 
the covering may be dispensed with under certain conditions. 

But in wind-tunnel model tests for the determination of the critical 
speed or the proof of flutter freedom, an externally geometrically and 
internally dynamically similar model is necessary. This problem calls 
for a model scale that is not too small and generally also a far-reaching 
geometrical similarity of the internal construction. When compressibility 
effects are involved, the ratio of the sonic velocities of the structural 
material and of the flowing medium must agree. In wind tunnels with air 
media this practically indicates aluminum as model material. 

It is difficult to meet all these requirements. The first dynamically 
similar models built in Germany were therefore still rather far from the 
dynamical similarity with full-scale design. Substantial improvements 
were attained by the tendency toward external and internal geometric simi- 
larity with the use of adhesive plastics such as Vinidur and adhesive thin 
aluminum sheets. They are available in the U.S.A. in finely gaged thick- 
ness for model designs. Even small forming rollers are utilized for such 
strips . Such models can actually be regarded as miniature copies of full- 
scale versions. Kinnaman (ref. 53) gives an insight into the model testing 
technique of the Boeing Co., which prefers a cheaper wood construction, 
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as was customary In Germany. Melzer (ref. 53a) gives an outline of the 
lavs of similarity and describes the technique of the Junkers Airplane Co. 

However , it doeB not always call for a complete model. Experiments 
with a half-span wing restrained at the root can be very informative, when, 
say, the effect of individual masses on the wing (engines, fuel tanks) on 
the critical speed is involved. Such simplified model tests have also 
been used as experimental check on corresponding three-dimensional flutter 
calculations, since the approximate assumptions, which must be introduced 
in order to simplify the theoretical calculation, can already be proved 
on simplified systems (Runyan-Sewall (ref. 54), Woolston-Runyan (ref. 55); 
Gayman (ref. 56 ), Tuovila, Baker and Regier (ref. 56 a) and Nelson and 
Rainey (ref. 56 b)). 

Since dynamically similar models of a complete airplane must not be 
too small, sufficiently great tunnels are required for their study, at 
airspeeds ranging within the permissible flying speeds of the full-scale 
design. This would be necessary if compressibility effects are to be 
accounted for. In favorable climate the expensive big tunnel can be 
replaced by a towing section in the open air, whereby a vehicle on rails, 
propelled by rockets, carries the model. Entire tail surfaces have been 
towed up to sonic velocity on such sleds, at the Edwards Air Force Base 
(U.S.A) and tested for flutter. A slight breeze, ground effect, and the 
inevitable vibrations are less disturbing in such tests than in steady 
profile measurements. 

Flutter studies have also been made on dynamically similar models in 
free flight (ref. 5). The models can be studied in diving flight or can 
be equipped with their own rocket drive and launched seaward in flight 
without return. The salient data are generally radioed to the ground 
where they are recorded. Such free-flying models have been used repeatedly 
in the U.S. for aeroelasticity problems as well as other important flight 
characteristics. Free-drop tests of models on which flutter begins 
directly above sonic velocity have been described by Dat, Destruynder, 
Loiseau, and Trubert (ref. 57) &rd compared with theoretical calculations. 


13. FLIGHT TESTS AT FUUL SCALE ' 


The flight safety of new airplane types must be proved by flight at 
maximum speed. Then it may happen that the critical speed is exceeded 
without initiating flutter, due to an accidentally absent outside impulse. 
On the other hand, sudden incipient flutter may lead to the destruction 
of the airplane before the pilot has time to observe the particulars of 
the event. Attempts have therefore been made to find some way by which 
a dangerous flutter possibility, that escaped theoretical or experimental 
detection, could still be spotted early enough in a flight test. 
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In Germany they mounted, for this purpose, oscillating masses on 
levers into the wing, which were excited with slowly increasing frequency, 
as in the static oscillation test. It was hoped that, at increasing 
oscillation amplitudes with increasing flying speed, an approach to the 
critical flutter oscillation could be identified at the right time. This 
is correct in principle. But in the execution of the flight tests it is 
very difficult and time-consuming to grade frequency and speed fine enough. 
The rise in amplitude can therefore be so fast that it is too late for a 
warning and the airplane crashes. However, favorable flutter modes whose 
amplitudes in a given speed range do not exceed a small, still nondangerous 
amount have been detected by this method. These were carefully recorded 
by flight tests, because they are suitable for checking theoretical methods 
of calculations . 

On airplanes with servocontrol, it is easier to impress a periodic 
booster force on the control surface which is then transformed into a 
many times greater periodic wing load distribution by aerodynamic trans- 
formation. It is not difficult to install an additional corresponding 
servocontrol for acceptance testing. This method has been used in the 
U.S.A. Pepping (ref. 58 ) made an analytical study of wing flutter by 
artificially excited control surface oscillations. The flutter stability 
can be determined from the initiated torque of the oscillat ing control 
surface per unit of angle of rotation and observed periodically; when it 
is zero, flutter condition is reached. This dangerous condition can be 
avoided by installing a hydraulic damper, or even more effectively, a 
feedback coupling with the wing rotation in the servocontrol. This raises 
the critical speed, and makes it possible to draw conclusions about the 
critical speed which the aircraft would have without stabilizing servo 
mechanism from the measurements of the amplitudes and phase angles of the 
excitation and the wing motion. When the rate of reaction of the servo 
mechanism is high enough, it can be utilized to compensate any potential 
flutter by appropriate control motions with appropriate feedback. Even 
gust stresses can be compensated this way. Admittedly, it may seem risky 
to rely on the functioning of a sensitive servo mechanism for suppressing 
flutter. But it can be used to good advantage, at least in high-speed 
tests of a new airplane type. 


14. PREVENTION OF FLUTTER BI DESIGN SPECIFICATIONS 


In spite of the multitude of measurements of the stress frequency 
distributions of airplanes in flight already available, it has seldom 
been possible to use the physical process of the flight stress as basis 
of the strength calculations of airplanes. The strength of modem air- 
craft in flight rests rather on very primitive dimension-analytical 
formulas and empirical values on conventional aircraft compiled in the 
design specifications for airplanes. 
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[Therefore, it seemed justifiable to follow a similar path in the 
practical treatment of flutter in order to be able to build airplanes 
with less flutter risk at tolerable expense. Ihe available theoretical 
and experimental aids are inadequate for a highly probable correct pre- 
diction of flutter freedom. 

Ihe convenient method developed in the early days of aviation, of 
piloting the airplane by free -swinging control surface hinged at the wing . 
trailing edge, actually promoted flutter. Ihe subsequent addition of 
auxiliary control surfaces multiplied the difficulty. Ihis fundamental 
defect of all airplanes up to now can be largely removed by placing the 
centroidal axis of the control surface and auxiliary control surface in 
or directly before their axis of rotation. Ihis way, the degrees of 
freedom "bending or translation of the wing" and "rotation of control 
surface" are uncoupled at least as regards the mass forces and Kelvin's 
impulses . The coupling of the degrees of freedom "wing rotation, " 

"control surface rotation, " and "auxiliary control surface rotation" is 
reduced also, but not eliminated. Complete uncoupling of the degrees of 
freedom is impossible for aerodynamic reasons. 

Ihe beneficial effect of this mass balance of control surfaces is 
voided when the control surface or auxiliary control surface with its 
control cables has static natural oscillations which are in resonance 
with those of the wing or with one another or approach it. Ihe phase 
angle between wing oscillation and control surface oscillation is then 
very slight, i.e., variable by very slight outside forces, so that it 
can assume the value most favorable for energy absorption from the air- 
stream and so lead to flutter. The second design rule therefore specifies 
the avoidance of such neighboring frequencies, the permissible clearance 
being defined by experiment and simplified flutter calculations. 

Ihe stepwise application of these construction rules for the design 
of new aircraft and their check on the finished * airplane by static oscil- 
lation tests lowered the flutter probability considerably in Germany during 
the 1933 to 1944 period, where the effect of the individual measures could 
be proved statistically. Table I gives the cases of flutter investigated 
In Germany according to the statistics of Kussner (ref. 59) and Schwarzmann 
(ref. 60), for the periods a = 1925-1933, b = 1954-1940, and 
c = 1941-1944. Every spontaneously occurring case of 'flutter is counted 
once, arbitrarily inducible flutter altogether counted only once. 

In these accidents more than 60 pilots and test engineers lost their 
lives, 4l of them due to auxiliary control surface flutter. Ihe surpris- 
ingly high proportion of the auxiliary control surfaces is due to the fact 
that these small but vital structural components and their control cables 
and supports were not always given the necessary care in the design and man- 
ufacture. By themselves positively moving, sufficiently stiff auxiliary 
control surfaces require no mass balance. But if they or their connections 
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axe too weak ox break, a new degree of freedom coupled by mass forces 
is produced. Hence, the requirement for mass balance of auxiliary con- 
trol surfaces in the German specifications. 

Looking back, it is apparent that a careful application of the rules: 

1. Displacement of axis of gravity in or directly before its axis of 
rotation 

2. Avoidance of frequencies approaching static oscillations on wings, 
control surfaces, and auxiliary control surfaces 

would have prevented most of these flutter incidents. Experience indi- 
cates further that the permissible departures from these rules, i.e., 
the manufacturing tolerances of production aircraft, must be very narrow. 
Changes that may seem slight, even a new coat of paint etc., have caused 
flutter on proved airplane types, more than once. 

For 39 of the observed cases of flutter calculations of the reduced 
frequency 



are available, where 2Z m is the mean wing chord of the wing of tail- 

surface portion under maximum flutter. Of particular interest are those 
values which are only rarely exceeded. Table II represents the highest, 
second highest, and third highest values. The highest values are of 
lesser interest, since the individual values may contain relatively great 
observation errors. The values = 2.4 was observed fairly accurately 

on a flying boat and represents a special case, because the aileron 
fluttered with two degrees of freedom. 

The physical basis for the empirical fact that certain maximum values 
of the reduced frequency are rarely exceeded must be looked for in the 
structural damping of the aircraft. The damping phase angle on conven- 
tional aircraft is, > 0.01. To be able to absorb the corresponding 

energy from the airstream requires certain minimum values of the wave 
length of the oscillatory motion even for optimum phase angle of degrees 
of freedom. On assuming that flutter frequency is always higher than 
the fundamental frequency of the static oscillation, flutter should 

be avoidable when 


v 



( 110 ) 


where o>j_ indicates one of the values of table II and Vq the maximum 
speed . 
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Originally, it was attempted in Germany to design airplanes according 
to equation (110) . Since the degrees of freedom "wing bending" and 
"aileron rotation" have been decoupled since 1933 by the mass balance 
recommended in the German design specifications, the fundamental wing 
rotation served as lowest flutter responsive static oscillation mode and 
a corresponding torsional wing stiffness was demanded. But this require- 
ment also proved itself as unnecessarily far-reaching for high-speed air- 
craft on which the two aforementioned design rules had been fulfilled, 
as numerous flutter calculations and model tests indicated. Auxiliary 
control surface flutter, which chiefly occurs with the degrees of freedom 
"control surface rotation" and "auxiliary control surface rotation" could 
not be prevented by wing stiffness no matter how great-. For these reasons, 
the overall frequency (110) was stricken from the German Specifications 1936 
(ref. 6l) . In its place a theoretical or experimental proof of flutter 
freedom in the permissible flight range was demanded . . This was still 
afflicted with many uncertainties at that time, but it exerted, at least, 
an educational effect on the airplane builder and clarified the influence 
of the most important flutter parameters. 

Through this flutter proof, which extended to the static limiting 
cases of wing deflection and reversal of control surface effect, the wing 
stiffness and control surface stiffness were given a lower limit. 

It was found, repeatedly, that the stiffness prevailing for reasons 
of strength itself was above this lower limit, hence, that the parti- 
cular airplane was safe against flutter and static instability when it 
satisfied the strength requirements and in addition satisfied the two 
aforementioned design regulations . Wittmeyer (ref. 62) advocated a 
corresponding airplane development while retaining certain basic design 
modes, which have proved themselves as flutter -safe in an extended series 
of developments. By prescribing certain dimensionless parameters and 
additional rules for the auxiliary control surfaces in the design of new 
types of aircraft, it should then be possible to attain flutter safety 
by way of compliance with the strength requirements. However, such design 
rules can be applicable only in a narrowly limited empirical range. When 
advancing into higher speed ranges and with the use of very thin airfoils, 
the matter becomes different. In section 9> there is given a third design 
rule which prevents skin buckling in any operating condition; de Vries 
(ref. 63 ) gives a critical summary of the design rules for the prevention 
of flutter. 1 

In the British Specifications and in those of the ICAO (international 
Civil Aviation Organization (ref. 6k)), overall stiffness requirements 
still take up a lot of space. It deals with recommendations for conven- 
tional aircraft, the compliance of which dispenses with a further flutter 
proof. These requirements call for minimum values of the mean torsional 
stiffness of wing, fuselage, tail unit, and control surfaces, which depend 
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upon the permissible maximum speed and on certain dimensions of these 
structural components as follows 


M > k 2 
0 “ 



°oJ 


£>£ 0.8 

c 0 


( 111 ) 


M is a moment, 0 the angle of rotation during this moment at a certain 
point, v-q the maximum gliding speed at 30° path slope, B and S are 
lengths; on the wing, for example 


B = 0.9b S « 


b = wing span from root to tip, the mean geometric chord. The 

dimensionless constant k for the separate components is specified. 

This formula (ill) is applied also to the flexural stiffness of the 
fuselage and the elevator overhang arm, with M and 0 defined 
accordingly . 

Equation (ill) follows from the condition of the static torsional 
stability of the -wing by forming the mean value of the integral. For 
this case, the Prandtl factor is correct. The constant k 2 iB then, 
however, proportional to the rearward position of the elastic axis^ 
behind the neutral axis at quarter wing chord, but this rearward posi- 
tion is materially dependent upon the respective structural design of 
the wing. 

So far as the overall stiffnesses according to equation (ill) ate 
to serve for the prevention of flutter, they are supportable only by 
equation (110) and table H, hence are disguised frequency requirements. 

(Cf. Collar, Broadbent, and Puttick (refs. 2 and 65 ).) But the static 
oscillation frequencies of an airplane are dependent in a rather 

complicated manner upon the size and distribution of masses and stiffnesses. 
To justify overall stiffness formulas like equation ( ill ), the constructive 
tolerance, conceded to the "conventional" airplanes, would have to be very 
narrow, i.e., they all would have to be practically geometrically similar 
and partially dynamically similar, according to development series inves- 
tigated by Wittmeyer (ref. 62). Or else the safety factor of the overall 
stiffness would have to be very great, since the actual differences of the 
airplanes of different manufacturers cause a correspondingly great varia- 
tion of the decisive oscillation frequencies v^. 

-'-The calculating methods hitherto were based on the replacement of the 
wing by an elastic beam with elastic axis. A plate-like wing has no 
elastic axis. 
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Shifting the gravity axis and the elastic axis Of the wing in the 
neutral axis, the degrees of freedom "wing bending" and "wing rotation" 
in the subsonic range can be decoupled, thus avoiding flutter. A mass- 
balanced control surface does not change this action very much as a rule. 
This ma86 balance of the wing has been used very rarely in the past. 

But it might be that future aircraft with very thin wings will be subject 
to a requirement regarding the position of the gravity axis of the wing, 
in order to lower the required torsional stiffness of the wing to a 
structurally tolerable level. 

To assure safety of aircraft against flutter in the future, the young 
aeronautical engineers must be made familiar with aeroelastic problems 
now, so as to enable them to take the correct safeguards while the air- 
craft is being designed. Formerly, airplanes were generally designed 
or even built without regard to flutter hazards, and trying to find a 
remedy proved then a too long and expensive undertaking. The flight 
performances often deteriorated so much that the particular type had to 
be withdrawn from the competition. A subsequent elevator mass balance, 
for example, cost 10 percent of the pay load and more on balance weight, 
aside from the greater air resistance of the externally mounted balance 
weights . 

As we enter the higher speed ranges, aeroelastic problems muBt 
receive particular attention. Empirical data and design rules, which 
are satisfactory or still acceptable for conventional aircraft, can no 
longer be relied upon, the physical process itself must be studied and 
mastered . 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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TABLE II.- MAXIMUM VALUES OF OBSERVED BEDUCED FREQUENCY 



°L 

“2 

*5 

Wing unit 

(2.4) 

1.29 

l.l4 

Horizontal tail surface 

1.1 

0.55 

0.51 

Vertical tail surface 

0.4- 

0.4-0 

0.58 
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